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Abstract. - Wo demonstrate direct evidence for ferromagnetic order at defect structures in 
highly oriented pyrolytic graphite with magnetic force microscopy at room temperature. Magnetic 
impurities have been excluded as the origin of the magnetic signal after careful analysis supporting 
an intrinsic magnetic behavior of carbon-based materials. The observed ferromagnetism has been 
attributed to originate from unpaired electron spins localized at grain boundaries. Scanning 
tumicling spectroscopy of grain boundaries showed intense localized states and enhanced charge 
density compared to bare graphite. 
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Introduction. — Graphite has been considered as a 
diamagnetic material for a long time. However, recent 
experiments have shown that ferromagnetic order is pos- 
sible in different carbon-based materials. Ferromagnetism 
with high Curie temperature, well above room tempera- 
ture, and very small saturation magnetization has been 
reported in various graphitic systems [1-6]. The role 
of different magnetic impurities on the measured ferro- 
magnetism has been studied in various samples of highly 
oriented graphite (HOPG), Kish graphite, and nature 
graphite [5]. The magnetization results showed no correla- 
tion with the magnetic impurity concentration [5] . Ferro- 
or ferrimagnetic ordering was demonstrated in proton- 
irradiated spots in highly oriented graphite [4]. Bulk fer- 
romagnetic graphite with a high defect concentration has 
been prepared via chemical route reaching the saturation 
magnetization 0.58 emu/g [3]. Apart from that, ferromag- 
netism has been observed in other carbon-based materials 
such as polymerized fuUcrcncs [7], carbon nanofoam [8], 
proton irradiated thin carbon films [9] and nitrogen and 
carbon ion implanted nanodiamond [10]. All these ob- 
servations suggest an inherent ferromagnetic behavior of 
carbon-based materials. 

Several theoretical investigations have been carried out 
to explain magnetism observed in these systems. The ori- 
gin of ferromagnetism was suggested to be attributed to 
the mixture of carbon atoms with sp^ and sp"^ bonds re- 
sulting in ferromagnetic interaction of spins separated by 
sp^ centers [11]. Another theoretical calculation suggested 
magnetism in sp^ bonded carbon nanostructures, which 



contains a negatively curved graphitic surface introduced 
via the presence of seven- or eight-membered rings [?]. 
In nanometer scale graphite, the electronic structure is 
strongly affected by the structure of the edges. Fujita 
and coworkers proposed that the tt electrons on a mono- 
hydrogenated zigzag edge might create a ferrimagnetic 
spin structure on the edge [12]. Recently, it has been 
shown in spin-polarized density functional theory calcu- 
lations that point defects in graphite such as vacancies 
and hydrogen-terminated vacancies are magnetic [13,14]. 
Three-dimensional network of single vacancies in graphite 
developed ferrimagnetic ordering up to 1 nm separation 
among the vacancies [30]. 

In this Letter, we report an experimental study of fer- 
romagnetic order in highly oriented pyrolytic graphite 
(HOPG) arising from defect structures. A ferromagnetic 
signal has been observed locally with magnetic force mi- 
croscopy (MFM) and in the bulk magnetization measure- 
ments using superconducting quantum interference device 
(SQUID). Observed ferromagnetism has been attributed 
to originate from itinerant electrons occupying narrow de- 
fect states of grain boundaries in the graphite crystal. 
A systematic study of the grain boundaries have been per- 
formed on the same samples with scanning tunneling mi- 
croscopy (STM) and spectroscopy (STS). 

Experimental. — Samples of HOPG of ZYH qual- 
ity were purchased from NT-MDT. The ZYH quality of 
HOPG with the mosaic spread 3.5° - 5° has been chosen 
because it provides a high population of step edges and 
grain boundaries on the graphite surface. HOPG sam- 
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pies were cleaved by an adhesive tape in air and trans- 
ferred into a scanning tunneling microscope (Omicron LT 
STM) working under ultra high vacuum (UHV) condition. 
The HOPG samples have been heated to 500° C in UHV 
before the STM experiments. STM measurements were 
performed at 78 K in the constant current mode with me- 
chanically formed Pt/Ir tips. The same samples have been 
subsequently studied by atomic force microscopy (AFM), 
magnetic force microscopy (MFM) and electrostatic force 
microscopy (EFM) in air using Dimension 3100 SPM from 
Veeco Instruments. PPP-MFMR cantilevers made by 
NanoSensors with a hard magnetic material Co-coating 
film have been used in the MFM tapping/lift mode. 

Results and discussion. — In figure 1, AFM, MFM 
and EFM images of the same area on the HOPG surface 
are shown. AFM topography picture in figure la displays 
a surface with a high population of step edges, surface 
distortions and defects. The MFM images in figure lb 
and figure Ic were taken on the same place as the AFM 
image with a lift scan height of 50 nm, where long-range 
van der Waals forces are negligible and magnetic forces 
prevail. Magnetic signal is measured on most of the line 
defects, however, a step edge marked as A in figure la does 
not show a magnetic signal in the MFM image. On the 
other hand, two lines in the MFM image in figure lb that 
are indicated as B and C do not show a noticeable height 
difference in the topography. The lines B and C are grain 
boundaries of HOPG. A detailed AFM and STM study 
can be found elsewhere [16,17]. 

In order to determine the character of the detected mag- 
netic signal, the MFM tip has been magnetized in two op- 
posite directions: aiming into (figure lb) and out of the 
graphite surface plane (figure Ic). Since the MFM signal 
represents the phase shift between the probe oscillation 
and the driving signal due to a magnetic force acting on 
the tip, the dependence of the phase shift on the force 
gradient can be expressed by a simple form [18] 



QdF^ 
k dz ' 



(1) 



where Q is quality factor and k is spring constant of the 
cantilever. Typical values of our MFM system give a min- 
imal detectable force gradient in the order of 100 /iN/m, 
Q — 200 and k — 2.8 N/m. For a true quantitative inter- 
pretation of MFM images it is necessary to have an exact 
knowledge of the geometry and magnetic properties of the 
tip and the substrate in order to express the force acting 
on the tip, which is difficult and has been achieved only in 
special cases [18]. Nevertheless, a qualitative analysis can 
be done according to expression 1, where a positive phase 
shift (bright contrast) represents a repulsive force between 
the tip and the sample, and a negative phase shift (dark 
contrast) manifests an attractive interaction relative to 
the background signal. Since the tip magnetized into the 
graphite surface plane has shown a bright contrast in fig- 
ure figure lb and out of plane magnetized tip produced a 




Fig. 1: (Color online) The same area on the HOPG surface 
imaged with AFM (a), MFM (b) and (c), and EFM (d). MFM 
tip has been magnetized into the graphite surface (b) and out 
of the graphite surface (c), respectively. Image parameters: 
san area 2x2 /im^, AFM z-range z — ^ nm, MFM z-range (b) 
$ = 2° and (c) $ = 1°, the MFM lift height /i = 50 nm, EFM 
z-range $ = 1°, the EFM lift height /i = 20 nm. 



dark phase contrast on the line defects in figure lb, the ori- 
entation of the net magnetic moment in the defects stayed 
in the same direction, pointing out of the graphite surface 
plane. This shows a clear indication of ferromagnetic or- 
der at the defect sites at room temperature. In the case of 
paramagnetic order, a bright contrast would be detected 
in both direction of the magnetization of the tip because 
the local magnetic moments would align with the mag- 
netic field of the tip leading to attractive interaction. The 
same result would be valid if electric force gradients were 
detected due to charge accumulation at the step edges. 
Therefore, the ferromagnetic order in the defects of the 
HOPG sample is the only plausible explanation for the 
detected MFM signal. 

However, not all the signal measured in the MFM 
showed to be sensitive to the reversal of the tip magneti- 
zation, in particular, areas with a different phase contrast. 
This is due to the metallic character of the magnetic coat- 
ing film of the MFM tip, which probes electrostatic forces 
as well. Therefore EFM has been measured on the same 
place with Pt coated Si tip with a lift scan height of 20 nm 
(see figure Id). A bright contrast is observed on the same 
places as in the MFM images. Similar observations of 
regions with a different potential has been measured in 
EFM and kelvin probe microscopy (KPM) on HOPG be- 
fore [31,32]. This non-uniform potential distribution has 
been found to be caused by the mechanical stress induced 
during sample cleaving [32]. Thereby the MFM measure- 
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Fig. 2: Out-of-plane (H \\ c) SQUID magnetization mea- 
surement on HOPG after substraction of the diamagnetic 
signal at 5 K (a) and 300 K (b). The diamagnetic back- 
ground signals were x ~ —3.2 x 10~* emu/g mT at 5 K and 
X = -1.0 X 10~* emu/g mT at 300 K. 

ments represent a superposition of magnetic and electro- 
static signal, which explains well the observed line shapes 
in figure Ic. 

The bulk magnetization of the HOPG samples have 
been analyzed with SQUID magnetometer at 5 K and 
300 K. In figure 2, out-of-plane magnetization loops of 
HOPG after substraction of linear diamagnetic back- 
ground signals are shown. They demonstrate small hys- 
teresis even at room temperature, which is a clear sign of 
ferromagnetic order. The out-of-plane saturation magne- 
tization reaches value 0.013 emu/g at 5 K. In-plane mag- 
netization loops have been measured on the same sample 
as well. They have shown ferromagnetic hysteresis loops 
comparable to SQUID measurements on HOPG reported 
by P. Esquinazi et al. [5]. The in-plane magnetization 
loops saturated at one order smaller values 2 x 10^"^ emu/g 
than in the out-of-plane configuration at 5 K. The coercive 
field and remnant magnetization have been found similar 
in both in-plane and out-of-plane magnetization measure- 
ments. In the work of P. Esquinazi at al. [5], the ferro- 
magnetic signals were measured up to temperature 500 K. 

The observed high temperature ferromagnetism in 
HOPG can have different possible origins. The first one 



is obviously ferromagnetism due to magnetic impurities. 
HOPG samples, as it has been studied previously [5], con- 
tain small fraction of magnetic elements. Therefore, we 
have analyzed the HOPG samples for impurity concen- 
tration by particle induced X-ray emission (PIXE) in the 
bulk material and by low energy ion scattering (LEIS) at 
the surface. As a main magnetic impurity in PIXE was 
found Fe with concentration around 20 fJ.g/g and Ti with 
concentration 4 /ig/g. Other magnetic and metallic impu- 
rities have been found below 1 /ig/g. The surface analysis 
by LEIS have not detected any magnetic elements indicat- 
ing that the concentration of all these elements is below 
100 ppm. The measured content of Fe impurities in HOPG 
is not sufficient to produce the ferromagnetic signal shown 
in figure 2. The amount of 1 /ig/g of Fe would contribute 
maximally 2.2 x 10""* emu/g to the magnetization and for 
Fe or Fe304 clusters, the magnetic signals would be even 
smaller [5]. 

Another possible source of the shown up ferromagnetic 
behavior are the defect structures in graphite. Line defects 
occur naturally in graphite as edges and grain boundaries. 
Graphite edges have been extensively studied both theo- 
retically [12,15] and experimentally [20-22]. There are two 
typical shapes for graphite edges: armchair and zigzag. 
Only zigzag edges are expected to give rise to the mag- 
netic ordering due to the existence of the edge state [12]. 
STM experimental results on step edges of graphite, how- 
ever, showed that zigzag edges are much smaller in length 
(« 2 nm) than those of armchair edges and less frequently 
observed [22]. At the same time, this would also sug- 
gest the occurrence of a one-dimensional ferromagnet of 
micrometer size at room temperature. We rather believe 
that step edges are created on HOPG at places where bulk 
grain boundaries cross the surface. During the cleavage of 
the graphite crystal these grain boundaries are the weak- 
est points of the graphite crystal. A step edge created by 
this way would have the same orientation and geometry 
as a grain boundary underneath. 

Grain boundaries in HOPG have been studied in great 
detail by AFM and STM [16,17]. Grain boundaries are 
inevitable defects in graphite because of polycrystalline 
character of the crystal. They are formed between two 
rotated grains during the crystal growth. They extend 
over step edges and form a continuous network all over 
the graphite surface. Grain boundaries show a small or 
no apparent height in AFM (see figure 1). On the other 
hand, grain boundaries exhibit a very distinct sign in 
STM, where they appear as a one dimensional superlattice 
with the height corrugation up to 1.5 nm. An example of 
STM and STS on a grain boundary with the periodicity 
13 = 1.4 nm is displayed in figure 3. Since grain bound- 
aries have almost no height in AFM and the corrugation 
of STM is given by convolution of the topography and the 
local density of states (DOS) of the substrate, hence grain 
boundaries must be endowed with an enhanced charge 
density compared to the bare graphite surface. Apparently 
this is seen in the STS of a grain boundary (figure 3b), 
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Fig. 3; (Color online) (Color online) (a) STM image of a grain 
boundary on HOPG showing a ID superlattice with periodicity 
D = 1.4 nm. The angle between two graphite grains is 18° and 
the angle between the grain boundary and graphite lattice is 
9°. Scanning parameters: 10 x 10 nm^, U = 0.6 V, It = 
0.4 nA. (b) STS on the grain boundary and on the bare graphite 
surface (tunneling resistance 0.9 GQ). The grain boundary 
shows localized states at -0.18 V and 0.36 V. 



which in addition shows two strong locaHzed states at - 
0.18 V and 0.36 V. These states are not observed on the 
bare graphite surface. Similarly, two localized states have 
been measured for different grain boundaries [17]. 

Defects in graphene break the electron-hole symmetry, 
which lead to creation of localized state at the Fermi en- 
ergy and to the phenomenon of self-doping, where charge 
is transferred to/from defects to the bulk [25,28]. The 
self-doping of defects is in accordance with our experimen- 
tal observations showing enhanced charge density at the 
grain boundaries. Since graphene systems have low elec- 
tron densities at the Fermi energy, electron-electron inter- 
actions play an important role as the recent experiments 
showed [25,27]. In the presence local repulsive electron- 
electron interaction the localized states will become polar- 
ized, leading to the formation of local moments [28] . This 
has been illustrated in DFT studies of point defects in 
graphite such as vacancies and hydrogen-terminated va- 
cancies. These defects revealed to be magnetic having 
a local magnetic moment larger than IfiB [13, 14]. In 
graphene, the indirect RKKY interaction mediated via 
valence electrons between these local moments has been 
found to be always ferromagnetic due to the semimetal- 
lic properties of graphene [28] . The presence of ferromag- 
netism has been studied theoretically in the phase diagram 
of pure and doped graphene [26] . It has been shown that 
ferromagnetic order can be stabilized at low doping when 
the exchange coupling is sufficiently large and as well at 
higher doping due to a disorder which further stabilizes the 
ferromagnetic phase. In graphite, however, because of fi- 
nite density of states at the Fermi energy there is expected 
to be competing ferro- and antiferromagnetic coupling be- 
tween the moments. This gives an additional oscillating 
term with the oscillation period determined by the Fermi 
momentum of electrons (holes) [28,29]. In the DFT study 
of an 3D array of single vacancies in graphite, different su- 
percells containing single vacancy have been studied [30]. 



Ferrimagnetic order has been supported up to the distance 
1 nm among the vacancies, while 5x5x1 supercell (1.23 
nm separated vacancies) did not show a net magnetic mo- 
ment in graphite [30]. In graphene, the 5x5 supercell 
exhibited still a net magnetic moment of l.72fiB [30]. 

In a similar way, grain boundary in one layer of graphite 
can be visualized as a one-dimensional line of equidis- 
tantly distributed defects, where the superlattice periodic- 
ity gives the distance between the defect sites. The defects 
in grain boundaries are not single vacancies, for which a 
trigonal symmetry would be expected to be observed in 
STM, but rather more complicated defects. If the spin 
polarized electron states were created in a grain bound- 
ary, the exchange splitting would be in order of 0.6 eV in 
our experiment (figure 3b). We have observed two split 
localized states around the Fermi energy at different grain 
boundaries with periodicities below 4 nm. Similarly like 
in the DFT study of an 3D array of single vacancies in 
graphite [30]. Another supporting evidence that the fer- 
romagnetism originates from grain boundaries is the fact 
that grain boundaries and step edges are the only defects 
observed in STM experiments on graphite surface. No 
point defects have been detected on the graphite surface 
with STM. Secondly important aspect is the two dimen- 
sional character of the grain boundaries, which explains 
most of the features from MFM and SQUID measure- 
ments. 

We assume that grain boundaries are propagating along 
the c-axis of the graphite crystal creating 2D plane of de- 
fects. The distance between the defects in grain bound- 
aries is determined by the superlattice periodicity in the 
plane and by the interlayer separation 0.33 nm along the 
c-axis. As it was described before, step edges can be the 
manifestation of the grain boundaries buried underneath 
them. The ferromagnetic signal would then come from 2D 
grain boundary planes formed through the bulk crystal. 
Moreover, an infinitely extended 2D magnetic plane with 
in-plane magnetization is stray-field-free and therefore it 
can exist in the single-domain state [23]. Accordingly, in- 
plane magnetized grain boundary plane should show a sin- 
gle magnetic domain, which supports the observation of 
only one magnetization direction in MFM measurements. 
Due to crossings among grain boundaries, minimum en- 
ergy configuration would lead to magnetization aiming in 
the c-axis of HOPG. 2D character of grain boundaries is 
also in accordance with the higher saturation magnetiza- 
tion measured along the c-axis than along the in-plane 
direction. It is expected that most of the grain bound- 
aries have small tilt towards the c-axis. Therefore, a larger 
magnetic field will be needed to align the local magnetic 
moments along the c-axis than in the graphene planes, 
where the magnetic axis stays in the 2D grain boundary 
plane. 

The electrons involved in the ferromagnetic behavior in 
graphite are sp electrons and therefore the well known 
theory of magnetism based upon the unfilled character of 
3d or 4/ electrons energy levels cannot be directly applied. 
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In the proton irradiated carbon studied by x-ray magnetic 
circular dichroism, it was demonstrated that the magnetic 
order originated only from the carbon 7r-electron system 
[9]. Similarly it has been found by Faccio et al. [30] that 
TT electrons are responsible for the magnetism in the 3D 
array of single vacanc;ies. 

Ferromagnetism with high Courie temperature Tq of 
itinerant sp electrons in narrow impurity (defect) bands 
has been recently studied in tlieory. This theory has been 
applied to CaBe but it is well applicable to graphite with 
narrow grain boundary states too. It has been argued that 
Stoner ferromagnetism with high Curie temperatures Tc 
can be expected for sp electron systems with narrow im- 
purity states [33] . In this theory, correlation effects do not 
reduce the effective interaction which enters the Stoner 
criterion in the same way as in a bulk band ferromagnets. 
Moreover, the spin wave excitations may not be effective 
if full spin alignment is maintained in lowering Tq- The 
value of Tc can be thus close to the value given by Stoner 
theory unlike for other bulk ferromagnetic metals. Ac- 
cording to the above analysis, we believe that the grain 
boundaries are playing a crucial role in the observed fer- 
romagnetism in HOPG. 

In conclusion, ferromagnetic signal has been observed 
in HOPG by magnetic force microscopy at room tempera- 
ture. The observed ferromagnetism has been attributed to 
originate from unpaired electron spins localized at defects 
sites of grain boundaries. Scanning tunneling microscopy 
(STM) and spectroscopy (STS) have revealed localized 
states and enhanced charge density of grain boundaries. 

* * * 

The authors are grateful to R. Lavrijsen and B. Koop- 
mans for SQUID measurements, P. H. A. Mutsaers for 
PIXE analysis and H. H. Brongersma for LEIS measure- 
ments. This research was supported by Nanoned. 

REFERENCES 

[I] KOPELEVICH Y., ESQUINAZI P., TORRES J. H. S. and 
MOEHLECKE S., J. Low. Temp. Phys., 119 (2000) 691. 

[2] COEY J. et al, Nature (London), 420 (2002) 156. 

[3] MOMBRU A. W. et al, Phys. Rev. B, 71 (2005) 100404. 

[4] EsQuiNAzi P., Spermann D., Hohne R., Setzer a., 

Han K.-H. and Butz T., Phys. Rev. Lett, 91 (2003) 

227201. 

[5] ESQUINAZI P. et al, Phys. Rev. B, 66 (2002) 024429. 

[6] ESQUINAZI P. and HoHNE R., J. Magn. Magn. Mater., 
290-291 (2005) 20. 

[7] Makarova T. L. et al. Nature, 413 (2001) 716; see how- 
over Retraction: Nature, 440 (2006) 707. 

[8] Rode A. V. et al, Phys. Rev. B, 70 (2004) 054407. 

[9] Ohldag H. et al, Phys. Rev. Lett, 98 (2007) 187204. 

[10] Talapatra S. et al, Phys. Rev. Lett, 95 (2005) 097201. 

[II] OvcHiNNiKOV A. A. and Spector V. N., Synth. Mat, 
27 (1988) 615. 



[12] FujiTA M., Wakabayashi K., Nakada K. and Kusak- 
ABE K., J. Phys. Soc. Jpn., 65 (1996) 1920. 

[13] Lehtinen p. O. et al, Phys. Rev. Lett, 93 (2004) 187202. 

[14] Yazyev O. V. AND Helm L., Phys. Rev. B, 75 (2007) 
125408. 

[15] Nakada K., Fujita M., Dresselhaus G. and Dressel- 
HAUS M. S., Phys. Rev. B, 54 (1996) 17954. 

[16] Cervenka J. and Flipse C. F. J., J. Phys.: Conf. Ser., 
61 (2007) 190. 

[17] Cervenka J. and Flipse C. F. J., (to be published). 

[18] Kebe T. and Carl A., J. Appl Phys., 95 (2004) 775. 

[19] Methling R.-P. et al, Eur. Phys. J. D, 16 (2001) 173. 

[20] Kobayasiii Y. et al, Phys. Rev. B, 71 (2005) 193406. 

[21] NiiMi Y. et al, Phys. Rev. B, 73 (2006) 085421. 

[22] KOBAYASHI Y. et al, Phys. Rev. B, 73 (2006) 125415. 

[23] Allenspacii R., J. Magn. Magn. Mat, 129 (1994) 160. 

[24] Wehling T. O. et al, Phys. Rev. B, 75 (2007) 125425. 

[25] Peres N. M. R., Guinea F. and Neto A. H. C, Phys. 
Rev. B, 73 (2006) 125411. 

[26] Peres N. M. R., Guinea F. and Neto A. H. C, Phys. 
Rev. B, 72 (2005) 174406. 

[27] GRliNElS A. et al, Phys. Rev. Lett., 100 (2008) 037601. 

[28] VozMEDiANO M. A. H., Lopez-Sancho M. p., Stauber 
T. and Guinea F., Phys. Rev. B, 72 (2005) 155121. 

[29] DuGAEV V. K., LiTViNOV V. I. and Barnas J., Phys. 
Rev. B, 74 (2006) 224438. 

[30] Faccio R. et al, Phys. Rev. B, 77 (2008) 035416. 

[31] Lu Y. et al, Phys. Rev. Lett, 97 (2006) 076805. 

[32] Marti'nez-Marti'n D. and Gomez-Herrero J., Cond- 
Mat, (2007) arXiv:0708.2994v2. 

[33] Edwards D. M. and Katsnelson M. L, J. Phys.: Con- 
dens. Matter, 18 (2006) 7209. 



p-5 



